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We predict a possible phonon softening instability in strongly correlated coupled semiconductor 
bilayer systems. By studying the plasmon-phonon coupling in coupled bilayer structures, we find that 
the renormalized acoustic phonon frequency may be softened at a finite wave vector due to many- 
body local field corrections, particularly in low density systems where correlation effects are strong. 
We discuss experimental possibilities to search for this predicted phonon softening phenomenon. 
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Much interest has recentlyEfcl focused on bilayer semi- 
conductor structures, where two quantum wells are put in 
close proximity with a high insulating barrier in-between 
to suppress interlayer tunneling. Such systems are of in- 
trinsic interest because competition between intra- and 
inter-layer Coulomb interaction (i.e., correlation effects), 
particularly at low carrier densities when Conktmb inter- 
action dominates (or in high magnetic fieldaju when the 
electronic kinetic energy is quenched), may lead to novel 
and exotic ground states, collective properties, and quan- 
tum phase transitions. Although much of the receat .work 
has concentrated on bilayer quantum Hall systems^ 13 the 
basic issue of the competition among electron kinetic en- 
ergy and intra- /inter-layer Coulomb correlations is quite 
generic, and interesting novel ground states and exotic 
collective properties are, in principle, possible even in the 
zero field situation provided the r s parameter defining 
the dimcnsionless ratio of Coulomb interaction strength 
to Fermi energy is high enough. In this Letter we pre- 
dict a particularly novel possibility, namely a finite wave 
vector acoustic phonon softening, driven by correlation 
effects in low density bilayer systems. We believe that 
our predicted phonon softening phenomenon may be re- 
alizable in the existing (electron or hole) bilayer systems 
provided the carrier density (and temperature) are low 
enough. 

Advances in nanofabrication and growth of semicon- 
ductor layered structures have recently made it possi- 
ble to fabricate very high quality (and very low density) 
GaAs-based high-mobility double well (bilayer) two di- 
mensional electron systems (2DES). The carrier density 
in these systems could be very low (large r s ), making 
them suitable for studying electronic correlation in low 
dimensions. As ±hc density decreases, unusual behav- 
ior is expectedcTEl due to the increasing importance of 
electron-electron interaction as the system goes from a 
weak-coupling (small r s ) to a strong-coupling (large r s ) 
regime. It is possible to have electron systems with r s 
as high as 15 - 20 and hole systems with r s as-iigh as 
25 - 30. One particularly interesting possibility^ is that 
an electronic charge density wave (CDW) instability may 
occur in these low density bilayer systems before the on- 
set of Wigner crystallization at still lower density. 

In bilayer structures there are two collective charge 



density excitations (plasmon modes^fl corresponding to 
the in-phase density fluctuation (optic plasmon, OP) and 
out-of-phase density fluctuation (acoustic plasmon, AP) 
between the layers. In the absence of any interlayer 
tunneling the AP mode has a linear dispersion in the 
2D wave vector because the long-range Coulomb interac- 
tion is suppressed by the screening from adjacent layers. 
Within random phase approximation (RPA) the acous- 
tic plasmon is shown to exist for small wave vector (long 
wavelength) at all densities. However, the linear disper- 
sion relation of the acoustic plasmon is sensitive to many- 
electron correlation effects even in the long wavelength 
regime. It has been shown that electron correlations may 
strongly renormalize the AP mode by pushing it com- 
pletely inside the single-particle electron-hole excitation 
(SPE) spectra, resulting in Landau damping of the made 
at much smaller wave vectors than predicted by RPAD'H. 
Thus, the AP mode may become overdamped due to cor- 
relation effects. Strong correlations can thus destroy the 
AP mode without smearing out of the boundary of the 
particle-hole continuum by defect or impurity scattering 
or thermal effectso. When this happens, there is a strong 
shift of single particle spectral weight toward zero energy 
at some finite wave vector due to the disappearance of 
the AP mode. The appearance of these spectral features 
with zero energy spectral weight at finite wave vector 
indicates that a state with a periodic modulation of den- 
sity costs relatively little energy to excite, and the system 
may therefore become unstable against- the formation of 
a periodic bilayer CDW ground stateta due to this AP 
plasmon overdamping mechanism. This CDW is differ- 
ent from the one associated with the 2kp (fcp being the 
Fermi wave vector) Peierls instability in one dimensional 
(ID) systemsSt2l, where the CDW instability occurs due 
to the electron-acoustic phonon interaction leading to a 
strongly divergent 2kp response function. The bilayer 
correlation induced CDW instability does not happen at 
2kp or any other Fermi surface related wave vector. 

In this Letter we study the electron-acoustic phonon 
interaction in low density bilayer systems. We investigate 
the plasmon (both OP and AP) mode coupling to the 
acoustic-phonons of the system and the consequent pos- 
sibility of a softening of acoustic phonon modes caused 
by the electron-phonon interaction in 2D bilayer systems. 
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The plasmon-phonon mode coupling phenomenon, which 
hybridizes the collective plasmon modes of the electron 
gas with the acoustic-phonon modes, gives rise to cou- 
pled plasmon-phonon modes. The mode coupling is par- 
ticularly strong at wave vectors around the plasmon- 
phonon resonance condition (u> p i ~ uj p h)- In 2D GaAs 
systems the resonant mode-coupling phenomenon nor- 
mally occurs at an extremely low density (n < 1.2 x 10 8 
cm~ 2 ) in the absence of correlation effects because the 
acoustic phonon energy is typically very low and one 
has to go to extremely low densities for the plasma fre- 
quency to be low enough to satisfy the w p i « u) p h condi- 
tion. Our important finding is that in the strongly corre- 
lated bilayer systems the plasmon-phonon mode coupling 
mechanism is sufficiently enhanced to give rise to strong 
(resonant) mode coupling effect at rather high densities 
which are already obtainable in currently available sam- 
ples (n < 1.0 x 10 10 cm~ 2 ). In addition, resonant mode- 
coupling occurs at large wave vectors (q > 2kp) mak- 
ing it potentially experimentally observable. The mode 
coupling effect depends strongly on electron density and 
layer separation since it originates from Coulomb correla- 
tions. We find that bilayer correlation effects in currently 
available samples could be strong enough that mode cou- 
pling may lead to a renormalized acoustic phonon mode 
which becomes soft at the wave vector around the reso- 
nance. The bilayer charge density instability associated 
with the overdamping of AP mode due to strong electron 
correlations gives rise to the drastic modification of the 
acoustic phonon dispersion leading to be softening. 

In this calculation we incorporate electron correlation 
effects in the local-field corrections which are deduced 
from quantum Monte Carlo numerical simulation data 
for the ground state of the 2DESLJ. We note that the 
specific correlation approximations used in our theory 
for the local field correction do not affect our qualita- 
tive conclusion about the existence of phonon softening 
in the low density bilayer systems, although the quan- 
titative details (precise density or wave vector at which 
softening occurs) depend on the approximation scheme. 
We take the bilayer system to be made of two identical 
symmetric quantum wells of width a separated by a dis- 
tance d. The Hamiltonian describing the bilayer electron- 
phonon system is given by H — Hq + H±, where Hq is 
the noninteracting electron and acoustic phonon kinetic 
energy, 

#o = E E E (*)4,J°k.i + E u Q h Q h Q> (*) 

k 1=1,2 Q 

and Hi is the electron-electron and the electron-phonon 
interaction, 

H i = E E v »'(<i) a l+ q ,i a l'- q ,i' a p'.i' a p,i 

k,k',q 1,1' =1,2 

+ E E A MQ)MqM- q j^Ab f Q + b- Q ), (2) 

k,Q 1=1,2 



where al dk,i are the l-th layer electron creation and de- 
struction operators with 2D wave vector k, Oq the phonon 
creation operator with 3D wave vector Q — (q,q z ), 
e(k) — k 2 /2m the kinetic energy of the electron (we 
take h = 1 throughout this paper), and u>q the bulk 
phonon frequency. Vw (q) is the Coulomb interaction be- 
tween l-th and l'-th layer whose matrix elements are given 
by Vw(q) = v(q)F w (q), where v^q) = 2ire 2 /(nq) with k 
as the background lattice dielectric constant and Fw (q) 
are the form factors describing the finite well width ef- 
fects. We consider only the lowest subband (for an in- 
finite square well potential) in each layer for simplicity. 
For two identical layers we have V\\ = V22 and V12 = V21. 
In Eq. (§), ft(q z ) = f dz^{z)* ^{z)e^\ where 0,(z) is 
the subband wave function in l-th layer. Mi(Q) is the 
coupling constant between the l-th layer electron and 
bulk GaAs phonons. In the jellium model Mi(Q) = 
i(2Tre 2 ujQ/Q 2 ) 1 / 2 . In GaAs/Al K Gai_ x As systems, elec- 
trons couple to acoustic phonons via deformation po- 
tential and piezoelectric couplings. We consider only 
the deformation potential coupling, which dominates in 
the present context. The electron-acoustic phonon cou- 
pling constant via the deformation-potential interaction 
is given by |M;(Q)| 2 = D 2 Q 2 / (2puiQ) where p is the mass 
density of the lattice, D the known deformation-potential 
coupling constant, and loq = coQ, where Co is the longi- 
tudinal bulk sound velocity, is the acoustic phonon fre- 
quency. 

Electrons interact among themselves through the di- 
rect Coulomb interaction and also through the dynamical 
virtual-acoustic-phonon exchange process. The effective 
interaction of an electron in l-th layer with an electron in 
l'-th layer is obtained in RPA by summing all the bubble 
diagrams 

Ww(q,u>) = ipw (q, u) + V>; m II(g, uj)W m i> (q, u>), (3) 

m 

where tpw (q, = V w (q) + Uw (<?, w), and the U w {q, w) is 
the phonon-mediated electron-electron interaction which 
is given by summing over the phonon wave vector in the 
^-direction (the 2D layers are in the x-y plane): 

Uu>{q^)=Y J M l (q,q z )Mt,{q,q z )\f l {q z )\ 2 D {Q,u), 

(4) 

where D°(Q,oj) = 2ojq/(oj 2 — oSq is the free-phonon 
propagator. We note that except for the form factor 
|/(g z )| 2 the integrand in Eq. (^) does not vary rapidly 
with q z . We therefore crudely approximate it with that 
factor multiplied by the integrand at q z = 0. (Explicit 
calculations show that this approximation is excellent.) 
Then we get 

Uw (q, w) = U ,° q \ [0(oj - Lu q )Ju> (qs) - 6(ojg - lo)I w (qs)} , 

(5) 
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where a(q,w) = ^/|1 -u 2 /u*\, and U (q) = D 2 q 2 /(2pc 2 ) 

for the deformation 

potential and Uo(q) = v(q) for the jellium model, respec- 



2 p -\z-z 



Jw{q) = 



tively, and I w {q) = J dzdz'\(j)i\ 2 \(t)i<\ 2 e 

The effective electron-electron interaction matrix, Eq 
(||), can be diagonalized through a canonical transfor- 
mation by introducing symmetric-antisymmetric (SAS) 

operators (ctk,± — [a>k,i ± Ofc,2]/v2) in the Hamiltonian. 
In the SAS representation the effective interaction then 
becomes decoupled by virtue of the symmetric nature of 
the bilayer system. The diagonalized elements of the ef- 
fective interaction are given by 



W±{q,u) 



l + ip±(q,uj)Il(q,w) 



v±(q) 



(6) 



where the plus (minus) label corresponds to the in- 
phase (out-of-phase) charge density modulation of the 
two layers. tp±(q,w) = V± + U±, where V±(q) — 
V n (q)[l - Gn(q)} ± V 12 (q)[l - G 12 (q)} is the Coulomb 
interaction matrix elements with the proper uxtealayer 
and interlayer local field corrections Gn^^uO; and 
U±(q,uj) — Un(q,uj) ± Ui2{q,ui) is the phonon mediated 
interaction matrix element; and H(q,u) the noninteract- 
ing 2D polarizability function. e t ^ t (q,ui) is the total di- 
electric function of the coupled electron-phonon system 
and is given by 



'(q^io) = l+V±(q)U(q,u) 



V±(q) + U ± (q,Lu)- 



(7) 



The dispersion relation of the coupled plasmon-acoustic 
phonon mode is given by the zeros of the dielectric func- 
tion e t ^ t (q,uj). We can rewrite Eq. (||) in terms of the 
screened Coulomb interaction and renormalized phonon 
interaction 



W ± (q,io) 



v±(q) 



D±{q,u). 



(8) 



The first term is the screened Coulomb interaction due 
to the collective motion of the electrons, and the second 
term represents the renormalized electron-phonon inter- 
action with the renormalized vertex M±/e± and phonon 
propagator D±. e±(q,u>) is the symmetric (antisym- 
metric) dielectric function of the electron system and is 
given by e±(q,u>) = 1 — V±(q)H(q,u>). The renormalized 
phonon propagator D±{q,uj) is 



D±(q,w) = 



2uj„ 



2uj q Ml(q,Lu)x±{q^Y 



(9) 



where X±(q> u> ) = ^-(q, u>)/e±(q,bj). The quasiparticle 
energies of the screened longitudinal vibrations are the 
poles of D{uj,q). The poles of the D{uj,q) are complex 
and hence we need an analytic continuation of the re- 
sponse function. However, it is a good approximation 
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FIG. 1. (a) The coupled plasmon-acoustic phonon modes 
and (b) the loss function (or the spectral weight) of the cou- 
pled mode, Im[l/e_ (q,u>)], as a function of energy for fixed 
wave vectors. Inset in (a) shows the coupled mode near the 
AP-acoustic phonon resonance. The area bounded by the 
light solid lines indicates the SPE Landau damping contin- 
uum. Arrows in (b) are the energies of the bare phonon mode 
at given wave vectors. 

to take the static limit (u> = 0) in the dielectric func- 
tion and to replace u> by the real frequency Q(q) in 
the imaginary part of the response function. Then we 
have the renormalized phonon modes Cl± and damping 



of the phonon 7± : 



2uj q MiRe X ±(q,uJ = 0), 



7±(9,w) = Mp.m X ±{q,n±(q))/[u q ~ 2M£Rex±(«,w - 
0)]. The width of the resonance in D±(q,uj) is given by 
27±(?>w). 

In Figs. 1 and 2 we show the results for the out-of- 
phase charge density modulation of the two layers in the 
bilayer system (the in-phase modulation, which is not 
shown, is similar to the single layer case). We use pa- 
rameters corresponding to electrons in GaAs quantum 
wellsEl: m = 0.067m e , k = 12.5, c = 5.14 x 10 5 cm/s, 
p = 5.3 g/cm 3 , and D = —16.0 eV for all our results. 
In Fig. 1 we show the coupled plasmon-acoustic phonon 
modes and the loss function, Im[l/e'°'(g, w)] which gives 
the mode spectral weight, as a function of energy for 
fixed wave vectors near the mode coupling region. In 
Fig. 1 we choose r s — 10 corresponding to the electron 
density n = 0.33 x 10 10 cm~ 2 , the well width a = lOOA, 
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and layer separation d — 300A. The area bounded by 
the light solid lines indicates the SPE Landau damping 
continuum. When correlation effects are turned off (i.e., 
G = 0) we have a well defined acoustic plasmon mode 
above the electron hole continuum and a well defined 
acoustic phonon mode with the bare phonon dispersion 
wph = coq. There is essentially no mode coupling without 
correlations because of the large difference in the plasmon 
and phonon energy scales. Local field effects give rise to a 
suppression of the AP mode above the SPE and its reap- 
pearance at large wave vectors beyond the single particle 
continuum in spite of it being overdamped (by Landau 
damping) at long wavelengths. This reappearance of the 
AP mode at large wave vectors, which is purely a lo- 
cal field effect, produces strong plasmon-phonon mode 
coupling. In Fig. 1(b) we show the spectral weight of 
the coupled modes. We find two peaks corresponding to 
the coupled modes, one below the bare acoustic phonon 
mode and the other above. Since the in-phase plasmon 
mode w+ (q) is essentially unaffected by correlation effects 
and has much higher energy than the phonon mode, the 
coupling of the in-phase plasmon mode to the acoustic 
phonon mode is negligible. (This is also the reason why 
the phonon softening phenomenon predicted in this Let- 
ter cannot occur in single layer 2DES.) 

In Fig. 2 we show the renormalized acoustic phonon 
mode dispersion and the damping for different layer sep- 
arations by finding the complex poles of the renormal- 
ized phonon propagator, Eq. (^J). The renormalized 
phonon mode is very different from the bare phonon 
mode due to correlation induced mode coupling. At the 
finite critical wave vector q c = 2.41fci? for d = 300A and 
q c = 1.75/cf for d — 690A we find phonon softening, 
where the acoustic phonon mode is renormalized to zero 
frequency u>(q c ) = 0. The critical wave vector depends on 
both density and layer separation, since it is a correlation 
effect. We emphasize that the bilayer phonon softening 
takes place due to the resonant coupling of GaAs acoustic 
phonons to the out-of-phase charge density modulation 
associated with the bilayer acoustic plasmons, which can 
become very strong at low enough carrier density. Strong 
acoustic phonon - acoustic plasmon coupling is the mech- 
anism causing the bilayer phonon softening. Our pre- 
dicted phonon softening thus necessarily requires both 
low carrier density (large r s ) and a bilayer system. 

We conclude by briefly discussing possible experimen- 
tal means to observe the predicted phonon softening. 
The most direct method would be to measure the bilaypt 
phonon dispersion using ballistic phonon spectroscopy^ 
as was done for single layer 2D systems some years ago. 
The other possibility is to look for phonon softening in- 
directly through bilayer drag measurementsEj where the 
magnitude of the drag should develop interesting anoma- 
lies in the soft mode regime since the drag is mediated 
partly by the phonon exchange mechanism. It may also 
be possible to .see the softening in bilayer transresistivity 
measurements^-]. In principle, it may be easier to see the 
softening phenomenon in hole bilayers where very large 
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FIG. 2. (a) The renormalized phonon modes £l p h(q) and 
(b) the damping of the phonon mode 7- (q, fi) as a function 
of wave vector for layer separations d = 300A (solid lines) and 
d = 690A (dashed lines) with r s = 10. We show the plasmon 
(AP) mode uj p i corresponding to d = 300A for comparison. 

r s values (r s « 20 — 30) can be achieved. 
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